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Abstract. We investigate pseudogap phenomena in the unitarity limit of a trapped superfluid 
Fermi gas. Including effect of strong pairing fluctuations within a T-matrix approximation, as 
well as effects of a harmonic trap within the local density approximation (LDA), we calculate 
the local superfluid density of states below the superfluid phase transition temperature T c . We 
show that the spatial region where single-particle excitations are dominated by the pseudogap 
may still exist even below T c , due to inhomogeneous pairing fluctuations caused by the trap 
potential. From the temperature dependence of the pseudogapped density of states, we identify 
the pseudogap regime of the unitarity Fermi gas with respect to the temperature and spatial 
position. We also show that the combined T-matrix theory with the LDA can quantitatively 
explain the local pressure which was recently observed in the unitarity limit of a 6 Li Fermi gas. 



1. Introduction 

Ultracold Fermi gases are very useful for the study of many-body physics in strongly 
interacting fermion systems. Indeed, using a tunable interaction associated with a Feshbach 
resonance [H [2J |3], we can study strong-coupling effects on single-particle properties, as well as 
thermodynamic properties[TJ HI E]j from the weak-coupling BCS (Bardeen-Cooper-Schrieffer) 
regime to the strong-coupling BEC (Bose-Einstein condensation) regime in a unified manner. 
In particular, the so-called pseudogap phenomenon, which is a typical phenomenon in strong- 
coupling fermion systems, has recently attracted much attention in cold Fermi gases. The 
photoemission spectra observed in a 40 K Fermi gas [6, [7j clearly deviates from those for a free 
Fermi gas in the BCS-BEC crossover region. As the origin of this, the importance of pairing 
fluctuations has been pointed out [9], [TU1 [TTJ [12l [13] . Since the pseudogap phenomenon has been 
also discussed in high-T c cupratesp3], the study of cold Fermi gases in the BCS-BEC crossover 
region is expected to be helpful for clarifying various anomalies observed in high-T c cuprates. 

In this paper, we investigate the local density of states (LDOS), as well as the local pressure, 
and effects of pseudogap effects associated with strong pairing fluctuations in a trapped unitarity 
Fermi superfluid. For this purpose, we employ a combined T-matrix theory with the local density 
approximation (LDA). We note that this strong-coupling theory has succeeded in quantitatively 



explaining the anomalous behavior of the photoemission spectra observed in a K Fermi 
gasJTQl ITT| 112] . Within this theoretical framework, we show that the spatial inhomogeneity 
due to the trap potential naturally leads to the coexistence of the spatial region where the 
ordinary BCS excitation gap can be seen and the region where the pseudogap still dominates 
over single-particle excitations in a gas cloud. We also apply this theory to the local pressure 
which has been recently observed in a 6 Li Fermi gas[5]. For simplicity, we set % = Jcb = 1 
throughout this paper. 



2. Combined T-matrix theory with local density approximation 

We consider a trapped superfluid Fermi gas with two atomic hyperfine states, described by 
pseudo spin a =f , j.. The model Hamiltonian in the absence of a trap is given by, in the Nambu 
representation [9] , 

^ = E*p^ 7 "3-At 1 ]* p -| £ Pj (q) Pj (-q). (1) 

P <JJ=1,2 

Here, ^> P = (c P ^,c_ p |) is the two-component Nambu field, where c Pa is the creation operator of 
a Fermi atom with the kinetic energy £ p = p 2 / (2m) — p, measured from the chemical potential 
fi. Tj(j = 1,2,3) are Pauli matrices acting on the particle-hole space, and the superfluid order 
parameter A is chosen so as to be proportional to the t\ component, pj = Y^p^ P + q T j^p(j = 
1,2) are the generalized density operators, describing amplitude (j = 1) and phase (j = 2) 
fluctuations in the Cooper channel. —U (< 0) is a pairing interaction which is assumed to be 
tunable. In this paper, we consider the case of unitarity limit, a^T 1 = 0, where the scattering 
length a s is related to the pairing interaction — U as = — U + J2p c ^-[15]- 

Effects of a trap is included within the LDA, which is conveniently achieved by replacing the 
Fermi chemical potential fi by the LDA expression fj,(r) = fj, — V{r), where V(r) = mco 2 r 2 /2 
is a harmonic trap. Various quantities then depend on the position r. For example, the 
LDA single-particle Green's function has the form G p (iw n ,r) ^ — Gp(iuj ri ,r) * — S p (iw n ,r), 
where G p (iw n ,r) _1 = iu n — C p (r)r3 + A(r)ri (where u) n is the fermion Matsubara frequency), 
£( r ) = p 2 /(2m) — /u(r), and A(r) is the LDA superfluid order parameter. S p (iw n ,r) is 
an LDA self-energy correction, involving effects of pairing fluctuations within the T-matrix 
approximation, which is given by S p (iw n ,r) = -T J2 q ,u n J2 s ,s'=± Fq S fan, r)r^ s G p+q (iu} n + 
ii/ n , r)r_ s /. Here, t± = t\ ± T2, and v n is the boson Matsubara frequency. T^ 5 (ii/ n ,r) = 
—U/[l + UH q s ' (iu n ,r)\ is the particle-particle scattering vertex, where H q s> (w n ,r) = 
T 2~2 q ,uj n T sG p+q / 2 (iu) n + iv n ,r)T s >G p _ q / 2 (iu n , r) is a pair correlation function. 

We self-consistently determine the superfluid order parameter A(r) and the chemical potential 
ji by solving the LDA gap equation, 



JUr) 2 + A(r) 2 

1 = 1 I tanh V^ j LL , ( 2) 

T 2,/£ P (0 2 + A(r)2 2T 



together with the equation for the number iV of Fermi atoms, 

POO 

N = 4irr 2 drn(r), (3) 

Jo 

where n(r) = 2TJ2p,u n G p (iui n , r)\ne lulnS is the particle density in the LDA. The superfluid 
phase transition temperature T c is determined as the temperature at which the gap equation 
([2]) with A(r) = is satisfied in the trap center r = 0|10j. We note that, below T c , the coupled 
equations ([2]) with Q give the situation that A(r) ^ when r > rg and A(r) = when r > tq, 
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Figure 1. Calculated local density of states p(uj,r) in the unitarity limit (T c = 0.312Tf, where 
Tp is the Fermi temperature), (a) T = T c . (b) T = 0.5T C . (c) T = 0.1T C . In each panel, we have 

offset the results, ro is obtained as (a) 0, (b) 0.64Rf> and 0.80-Rf, where Rp = \J 2ep / muQ is 
the Thomas- Fermi radius (where £f is the Fermi energy). Panel (d) is the phase diagram in the 
unitarity limit. "SF" is the region where the superfluid gap can be seen in LDOS. "PG" is the 
region dominated by pseudogap. "NF" is the region where LDOS is close to that for a normal 
Fermi gas. The solid line in panel (d) shows r = ro(T). 



where ro becomes large with decreasing the temperature. However, this is an artifact of the 
LDA, because A(r) should be finite everywhere in the gas below T c . Thus, we should regard ro 
as a characteristic radius, in the inside of which |A(r)| is large to some extent. 

Once A(r) and /j are determined, LDOS is calculated from p(uj,r) = —^J2p^ m G p (iu n — > 
u + iS, r) |ii- For the local pressure -P(r), using the relation dP{r) = n(r)dfi(r)^\W\ for a given 
temperature, we obtain 

= l9V{r) r r f dr f n r r fy ,q 
r or Joe 

3. Local Density of States and Local Pressure below T c 

Figure 0Ja) shows the LDOS in unitarity limit at T c . Although A(r) vanishes at T c , we see a 
large dip structure around u = in the trap center, originating from strong pairing fluctuations. 
However, such a dip structure is absent around the edge of the gas. Below T c , the pseudogap 
in the trap center is suppressed by the presence of finite A(r). Instead, the superfluid gap with 
the coherence peaks at the excitation edges appears (See Figs(TJb) an d ( c ))- 111 panel (b), the 

tanh ep(ro) 

pseudogap is found to appear around r = ro. By definition, 1 = UJ2 P 2 g (ro) * s satisfied at 
r = ro, which is the same form as the T c -equation when ro is replaced by r = 0. Thus, as in the 
case of the trap center at T c , strong pairing fluctuations induce the pseudogap around r = ro. 

When we define the pseudogap regime as the region where the dip structure appears in 
LDOS, we obtain the phase diagram shown in Fig. GD^d). As expected, the pseudogap regime 
(PG) exists around the "ro-line." We briefly note that the pseudogap regime exists even below 
this line where A(r) is finite. Below this pseudogap regime (denoted by "SF" in FigfTJd)), the 
superfluid gap accompanied by the coherence peaks appear in the LDOS. On the other hand, 
above the pseudogap regime, neither the pseudogap nor the superfluid gap appears in LDOS, 
which is denoted by "NF" in panel (d). Figure [ljcl) indicates that the SF region and PG region 
coexists below T c in a trapped superfluid Fermi gas. 

Strong-coupling effects can be also seen in the local pressure -P(r), as shown in Fig. [2J That 
is, -P(r) is remarkably enhanced from that for a free Fermi gas -Po( r ) (where the same chemical 
potential /i that used in -P(r) is also used.) In this figure, P(r) is found to be almost temperature 
independent, reflecting the university of the unitarity Fermi gas. As shown in this figure, our 
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Figure 2. Calculated local pressure P(r) in the unitarity limit. Solid circles are experimental 
results observed by ENS group [5j. P(r) is normalized by the local pressure Pq(t) in the case of 
a trapped free Fermi gas with the same chemical potential \x. 

results agree well with the recent experiment (solid circles) done by ENS group [5]. 
4. Summary 

To summarize, we have discussed inhomogeneous pseudogap phenomena in the unitarity limit of 
a trapped superfluid Fermi gas. Within the framework of the combined T-matrix theory with the 
LDA, we determined the spatial and temperature region in the gas cloud where the pseudogap 
appears. We also showed how the local pressure is affected by strong pairing interaction. The 
quantitative agreement of our results with the recent experiment by ENS group indicates the 
validity of the strong-coupling theory used in this paper for the superfluid Fermi gas in the 
unitarity regime. Our results would be helpful for the study of strong-coupling effects on trapped 
Fermi superfluids. 
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